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Abstract. We present a partial upgrade of the Monte Carlo event generator TAUOLA with the two 
and three hadron decay modes using the theoretical models based on Resonance Chiral Theory. 
These modes account for 88% of total hadronic width of the tau meson. First results of the model 
parameters have been obtained using BaBar data for 3n mode. 
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1. INTRODUCTION 

Tau lepton is the only lepton that decays into hadrons. From the perspective of low 
energies, as explored at B-factories, BaBar [1] and Belle [2] experiments, hadronic decay 
modes of tau lepton provide an excellent laboratory for modeling hadronic interactions 
at the energy scale of about 1-2 GeV, where neither perturbative QCD methods nor chiral 
lagrangians are expected to work with a good precision [3]. From the perspective of high 
energy experiments such as those at LHC, good understanding of tau leptons properties 
contributes important ingredients of new physics signatures. With the discovery of a new 
particle around the mass of 125-126GeV [4], tau decays are an important ingredient for 
determining if this is the Standard Model Higgs. This is especially pertinent since CMS 
has reported a deficit in the number of fermion decays from the new particle relative to 
the Standard Model Higgs Prediction. 

The first version of the program TAUOLA [5] was written in the 90's. TAUOLA 
is Monte Carlo event generator which simulates tau decays for both the leptonic and 
hadronic decays modes. The hadronic currents implemented in TAUOLA are based 
on resonance dominance model (RDM) [6] 2 . In the framework of RDM the hadronic 



Speaker 

2 Based on [5] Cleo and Aleph collaborations developed their own versions, non published, which differ 
from one another by the intermediate resonance states because of different detector sensitivity. 
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current for a three-pseudoscalar decay mode is a sum of weighted products of the Breit- 
Wigner functions. It is demonstrated in [7] that this approach is able to reproduce only 
LO #PT result. Moreover, the model was not sufficient to describe the Cleo KKn data 
[8]. This resulted in the Cleo collaboration reshaping the model by introducing two 
ad-hoc parameters that spoilt the QCD normalization for Weiss-Zumino contribution. 
The parameters are obtained fitting to data. However, before making conclusion that 
the Wess-Zumino anomaly normalization is spoilt it should be checked whether on 
oversimplified theoretical approximation, like RDM, was applied. 

The alternative approach based on the Resonance Chiral Theory (R^T) [9] was 
proposed in [7]. The computations done within R^T are able to reproduce the low- 
energy limit of ^PT at least up to NLO and demonstrate the right falloff in the high 
energy region [10]. The hadronic currents for the two and three meson modes (nn, Kn, 
KK, nnn and KKn) calculated in the framework of RjT were implemented in TAUOLA 
[11]. The tar-ball version of the code is available at web-page [12]. 



2. RESULTS FOR TWO AND THREE MESON FINAL STATES 

It is of utmost importance to implement the theory in a form that is as useful for 
general applications as possible. Therefore the corresponding hadronic currents (7^) in 
the upgraded version of TAUOLA have been written in the most general form compatible 
with Lorentz invariance. For x decay channels with two mesons (/zi(/?i) and h 2 (p 2 )) 

j^=N[{p l -p 2 -— (pi+P2)rF v ( S )+^(( Pl +p 2 rF s ( S )], a) 

s s 

where A12 = m \ — m\, s = (pi+ P2) 1 ■ For the final state of three pseudoscalars, with 
momenta p\,pi and p 3 , Lorentz invariance determines the decomposition 

= N {T^[c l (p2-p3) v F 1 +c 2 (p3-pi) v F2 + c 3 (p 1 -p 2 ) v F 3 ] 

+ C 4^ F 4 ~ 4^ c 5^vpoPlP2P3 F 5}, (2) 

where: T^ v = g^ v — q^q v /q 2 denotes the transverse projector, q^ = (p\ + p 2 + Pi)^ is 
the momentum of the hadronic system, F is the pion decay constant in the chiral limit. 

The decomposition (1) and (2) is the most general one, model-dependence is included 
in the hadronic form factors {Fy, F$ as well as F\, i = 1...5). The hadronic form factors 
calculated within R^T can be written as 

Fl = F X +F R +F RR (3) 
where Ff is the chiral contribution, Ff is the one resonance contribution and 

F RR 

is the double-resonance part. The explict form of the functions F{ for 3k and KKn 
modes can be found in [11], Section 2, as well as in [7, 13]. The theoretical results 
for the hadronic currents were obtained in the isospin limit (m K = (2m n + +m K o)/3, 
nix = (rn K + + m K o)/2), except for the two pion and two kaons modes. In the phase 
space generation, the differences between neutral and charged pion and kaon masses is 
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TABLE 1. The x decay partial widths. The PDG value [14] (2nd collumn) is compared with numer- 
ical results from TAUOLA with the RjT currents: in the isospin limit for pseudoscalar masses (3rd 
collumn), using physical masses (4th collumn). 



Channel 


Width, [GeV] 




PDG 


Equal masses 


Phase space 
with masses 


n~n° 


(5.778 ±0.35%) 


• 10 


-13 


(5.2283 ±0.005%) • 


io- 


-13 


(5.2441 ±0.005%) 


io- 


-13 


K-7t° 


(9.72 ±3.5%) 


io- 


15 


(8.3981 ±0.005%) • 


io- 


15 


(8.5810±0.005%) 


io- 


-15 


7t-K° 


(1.9 ±5% )■ 


io- 


14 


(1.6798 ±0.006%) • 


io- 


-14 


(1.65 12 ±0.006%) 


io- 


-14 


K-K° 


(3.60 ±10% ) 


io- 


15 


(2.6502±0.007%)- 


io- 


15 


(2.6502 ±0.008%) 


io- 


-15 


%~%~% + 


(2.11 ±0.8% ) 


io- 


-13 


(2.1013 ±0.016%) • 


io- 


-13 


(2.0800 ±0.017%) 


io- 


-13 


n°n°7t~ 


(2.10 ±1.2% ) 


io- 


-13 


(2.1013 ±0.016%) • 


io- 


-13 


(2.1256 ±0.017%) 


io- 


-13 


K-n-K+ 


(3.17 ±4% )• 


io- 


15 


(3.7379 ±0.024%) • 


io- 


15 


(3.8460 ±0.024%) 


io- 


-15 


K°K-K° 


(3.9 ±24% )• 


io- 


15 


(3.7385 ±0.024%) • 


io- 


15 


(3.5917 ±0.024%) 


io- 


-15 


K-7t Q K° 


(3.60 ±12.6% ) 


• 10 


-15 


(2.7367 ±0.025%) • 


io- 


15 


(2.7711 ±0.024%) 


io- 


-15 



taken into account, i.e. physical values are chosen. This has to be done to obtain proper 
kinematic configurations. The model parameters, more specifically the masses of the 
resonances and the coupling constants, were fitted to Aleph data, requiring correct high- 
energy behaviour of the related form factors, see Appendix C in [1 1]. 

To check stability of multidimensional integration in TAUOLA the MC results were 
compared with the semi-analytical ones (Gauss integration of analytical results was 
used). The difference between MC prediction and semi-analytical results for the partial 
decay width is less than 0.02% for all channels. Both differential spectra and numerical 
tests are collected at [12]. 

In Table 1, the partial decay width values from PDG [14] are compared with our 
results obtained with isospin-averaged pseudoscalar masses and with the physical ones. 
Comparison of the last two columns illustrates numerical effect of physical masses. The 
difference between the MC result and PDG one is 2% — 24% depending on the channel. 
As expected, that agreement is not good because only minimal attempts on adjusting to 
the model parameters have been applied for the comparison with BaBar and Belle data. 
The next section presents the first step toward this direction. 

3. FIT FOR THREE PION MODE TO BABAR DATA 

The main problem with upgrading the MC simulation of T decays is a lack of the 
published spectrum. Currently, only the two pion modes [2] and three pseudoscalar 
modes [15] are published. 

The results from a preliminary fit to the three pion mode (k + k~k~) [15] can be seen 
in Fig. 1 and Table 2. Both the three particle and the n + n~ invariant mass distributions 
have been considered. Disagreement about 12% level is visible in the low energy region 
of two particle invariant mass distribution whereas for the three invariant mass spectrum 
the difference between MC and data is less than 7%. However, the R^T parametrization 
is in better agreement with BaBar data than CLEO one, see Fig. 1. 

The fit was done taking into account only P-wave contribution of two pion system. 
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TABLE 2. Numerical values of the RjT parameters fitted to BaBar data for three pion mode 





M p , 


r p' 


M ai 


F 


Fv 


F A 


/V 


Min. 


1.44 


0.32 


1.00 


0.0920 


0.12 


0.1 


-0.36 


Max 


1.48 


0.39 


1.24 


0.0924 


0.24 


0.2 


-0.18 


Default 


1.453 


0.4 


1.12 


0.0924 


0.18 


0.149 


-0.25 


Fit* 


1.4302 


0.376061 


1.21706 


0.092318 


0.121938 


0.11291 


-0.208811 



* to[l5],X 2 /ndf = 2262.12/132 
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FIGURE 1. Invariant n + n~n~ (up) and n + n~ (down) mass distributions: the plots on the left-hand side 
correspond to the differential decay distribution, the ones on the right-hand side to plot ratios between MC 
and BaBar data. Lighter grey histograms are for RjT parametrization, darker grey is for CLEO one. 



As suggesting in [16] the discrepance in the low mass region could be described using 
a contribution from the scalar particle, S-wave contribution. We expect that inclusion 
of S-wave contribution [10] will improve the value of X 2 - The values in the fifth row 
of Table 2 are only the preliminary results. They will not necessarily correspond to the 
minimum of X 1 jndf of the final fit. Work is in progress. 
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4. CONCLUSION 



The theoretical results for the hadronic currents of two and three pseudoscalar modes, 
namely, %%, Kn, KK, nnn and KKn, in the framework of R^T have been implemented 
in TAUOLA. These modes, together with the one-meson decay modes, represent more 
than 88% of the hadronic width of tau. R^T is a more controlled QCD-based model 
than the usual used Breit-Wigner parametrizations. However, before making conclusion 
about validity of the model the theoretical results have to be confronted with the exper- 
imental data, which requires fit of the model parameters. Now that the technical work 
on current installation is complete, the work on fits is in progress in collaboration with 
theoreticians and experimentalists. 

At LHC tau decays are only used for identification and are not used to study their dy- 
namic. However, the dynamics of tau decays are important for both modeling the decays 
and -therefore the reconstruction and identification- and for measuring the polarization 
of tau decays. Therefore, an upgrade to the TAUOLA based on the BaBar and Belle 
results on tau decays is urgently needed for LHC. 

ACKNOWLEDGMENTS 

This research was supported by a Marie Curie Intra European Fellowship within the 
7th European Community Framework Programme (O.S.) and by Alexander von Hum- 
boldt Foundation (I.N.), by the Spanish Consolider Ingenio 2010 Programme CPAN 
(CSD2007-00042) and by MEC (Spain) under Grants FPA2007-60323, FPA201 1-23778 
(O.S. and PR.) and FPA201 1-25948 (PR.) and in part by the funds of Polish National 
Science Centre under decision DEC-2011/03/B/ST2/00107 (TP, Z.W.). 

REFERENCES 

1. B. Aubert et al., Phys.Rev.Lett. 100 (2008) 01 1801; J.P. Lees et al., arXiv: 1 109.1527 [hep-ex]. 

2. M. Fujikawa et al., Phys. Rev. D78 (2008) 072006; H. Hayashii, M. Fujikawa, Nucl.Phys.Proc.Suppl. 
198 (2010) 157. 

3. E. Braaten and C.-S. Li, Phys.Rev. D42 (1990) 3888. 

4. G. Aad et al., arXiv: 1207.7214 [hep-ex]; S. Chatrchyan, arXiv: 1207.7235 [hep-ex]. 

5. S. Jadach, Z. Was, R. Decker and J. H. Kiihn, Comput. Phys. Commun. 76 (1993) 361. 

6. J.H. Kiihn, E. Mirkes, Z. Phys. C56 (1992) 661; J.H. Kiihn, A. Santamarfa, Z. Phys. C48 (1990) 445. 

7. D. G. Dumm, A. Pich and J. Portoles, Phys. Rev. D69 (2004) 073002. 

8. T.E. Coan et al., Phys.Rev.Lett. 92 (2004) 232001 . 

9. G. Ecker, J. Gasser, H. Leutwyler, A. Pich, E. de Rafael, Phys. Lett. B223 (1989) 425; Nucl. Phys. 
8321(1989)311. 

10. P. Roig, I. Nugent, T. Przedzinski, O. Shekhovtsova, Z. Was, in these Proceedings. 

1 1. O. Shekhovtsova, T. Przedzinski, P. Roig and Z. Was, arXiv: 1203.3955 [hep-ph]. 

12. T. Przedzinski, O. Shekhovtsova and Z. Was, http://annapurna.ijj. edu.pl/^wasm/RChL/RChL.htm. 

13. D.G Dumm, P. Roig, A. Pich, J. Portoles, Phys. Lett. B685 (2010) 158; Phys. Rev. D81 (2010) 
034031. 

14. K. Nakamura, J. Phys. G37 (2010) 075021. 

15. I.M. Nugent, SLAC-R-936. 

16. E.I. Shibata et al., arXiv: hep-ex/0210039. 



MC generator TAUOLA: implementation of Resonance Chiral Theory for two and three meson modes. Comparison with experimer 



